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ABSTRACT 



Smart antennas possess the capability of suppressing jamming signal, so they can improve signal to 
interference plus noise ratio (SINR). Array processing utilizes information regarding locations of signal to aid 
in interference suppression and signal enhancement and is considered promising technology for anti -jamming. 
In this paper we studied three beam forming algorithms, Least Mean Square (LMS) algorithm, Optimized - 
LMS algorithm and Recursive Least Squares (RLS) algorithm. Simulation results are presented to compare the 
ability of these three algorithms to form beam in the direction of desired signal and place null in the direction 
of interference signal. Dependency of these algorithms on SNR and SIR is also analyzed. It has been found that 
RLS algorithm is best suited for anti -jamming applications. 
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1. INTRODUCTION 

Potential jamming in military and critical civilian 
applications has been a major concern for system 
designers. And usual filtering techniques are not 
helpful as the jamming signal and desired signal are 
of same frequency. Various methods have been 
adopted to avoid jamming; including frequency 
hopping but it requires excessive bandwidth. Spatial 
filtering can solve the problem [10] over head 
without the need of additional bandwidth as signals 
are filtered on basis of their direction of arrival. 



domain algorithms [6], and recently CSLMS 
algorithm [9]. Optimized-LMS algorithm modifies 
the conventional LMS algorithm with optimized step 
size and is studied in detail in this paper. The 
computational complexity of Optimized-LMS is also 
O (M). RLS algorithm usually converges with order 
of magnitude faster than LMS algorithm but the price 
paid is added complexity. Several variants of RLS 
algorithm are also proposed one of which is GVFF 
(Gradient based Variable Forgetting Factor)-RLS) [5]. 
The complexity of RLS algorithm is O (M2). 



Non-blind algorithms as discussed in this paper 
require the information of desired signal but blind 
algorithms such as Constant Modulus Algorithm 
(CMA) and MUSIC algorithm can estimate the 
Direction of Arrival (DOA) of the source signal, and 
then this direction information can be utilized in 
non-blind beam forming algorithms to form beam in 
the estimated direction. LMS algorithm is known for 
its simplicity and robustness. The computation 
complexity of LMS algorithm is O (M) . While it lacks 
in convergence speed several modifications to the 
algorithm are proposed including Optimized-LMS 
[2] Variable Step Size LMS (VSS-LMS) algorithms [6, 
7, 1, 3], variable-length LMS algorithm [7], transform 




Figure 1. Least Mean Square (LMS) Algorithm 
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In Fig. 1 the outputs of the individual sensors are 
linearly combined after being scaled with 
corresponding weights optimizing the antenna array 
to have maximum gain in the direction of desired 
signal and nulls in the direction of interferers. 

For beam former the output at any time k, y(k) is 

given by a linear combination of the data at M 
antennas, with u(k) being the input vector and w(k) 
being the weight vector 

y(k)=W H (k)u(k) 

M-l 

where w(k) = ^ 



k=0 



M-l 



and U (£) = £ 



u, 



k=0 



The signal received at time index k is 



+ u,(k- l)w, (k) + n(k) 



r(k) = u l (k — l)w l (m) + . . . 
= Y J u j (k- j)Wj +n(k) 

7=1 

= u T (k)W(k) + n(k) 



The output y(k) of the adaptive filter is expressed as 



y(k) -d x {k- l)w 1 (&) + ... 
^djik-fiwjik) 



.., + di (k -l)Wj (m) 



7=1 



= D T (k)w(k) 

In practice, the adaptive filter can only adjust w(k) 
such that y(k) closely approximates desired signal 
over time. Therefore, the instantaneous estimated 
error signal needed to up- date the weights of the 
adaptive filter is 

j(k) = p(k)e T (k)e(k) 
e(k) = r(k) - y(k) 

= r(k) - D T (k)w(k) 



This priori error signal, e(k) is used to minimize the 
estimator error by adaptive updating of filter 
weights. 



2.1 Least Mean Squares (LMS) Algorithm 

The LMS algorithm is based on the stochastic 
gradient and is given by [12] 

e(k) = u T (k)w(k) + zik) -u r (k)h(k) 
w(k + 1) = w(k) + r/u(k)e(k) 

where 77 is step size, u(k) is the transmitted diagonal 

matrix at sampling time k, W(k) is the adaptive filter 
coefficient, and e(k) is the estimation error. The filter 
coefficients are updated using an estimate of the cost 
function gradient, [pi(k)e(k) ■ In all practical 

applications, the signals involved might be corrupted 
by noise. When the noise is present in the received 
sequence, interference will also in the coefficients 
adaption process through the term \\\{k)e{k) As a 

result, where the distribution of the noise is highly 
impulsive, the LMS scheme might have low 
convergence and lower steady state MSE 
performance. The step size parameter, 77 determines 

the convergence rate of the algorithm and higher 
value provides faster convergence. However, if h 
exceeds certain bound then the algorithm will 
diverge. As the bound on h is not known a priori and 
is dependent on the various statistics 

2.2 Normalized LMS (NLMS) Algorithm 

The main problem of the LMS CE algorithm is that it 
is sensitive to the scaling of its input signals. This 
makes it very hard to choose h that guarantees 
stability of the algorithm. The NLMS is a variant of 
the LMS algorithm that solves this problem by 
normalizing with the power of the input signal. The 
NLMS algorithm can be summarized as [11]: 

w(k + 1) = w(k) + 7je(k) \ T (k)u(k) ^u(k) 

when a constant scalar step size is employed in the 
LMS are NLMS algorithm, there is a trade off among 
the steady state error-convergence towards the true 
channel coefficients, which avoids a fast convergence 
when the step size is preferred to be small for small 
output estimation error. In order to guarantee the 
algorithm to be convergent, the range of step size is 
specified but the choice of optimal learning step size 
has not been appropriately addressed. In order to 
deal with these troubles, one key idea is to exploit 
varying step size during adaptation. 
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2.3 Recursive Least Squares (RLS) Algorithm 

To combat the channel dynamics, the RLS based CE 
algorithm is frequently used for rapid convergence 
and improved MSE performance. The standard RLS 
algorithm is 

e(k) = u T (k)w(k) + z(k) - u T (k)w(k) 

R(k) = r l P(k - 1) - r l R(k)u T (k)P(k - 1) 
w(k + 1) = w(k) + xx(k)e{k)R{k) 

Where X is the exponential forgetting factor with 0 < 
X<1. The smaller value of 1 leads to faster 
convergence rate as well as larger fluctuations in the 
weight signal after the initial convergence. On the 
other hand, too small 1 value makes this algorithm 
unstable. Subsequently, it requires best possible 
forgetting factor such that the estimator error is 
decreased. 

Although a lot of modified CE algorithm has been 
studied on employing adaptive forgetting factor and 
parallel forgetting factor, the CE performance is 
severely degraded in highly dynamic fading channel 
even when the forgetting factor is well optimized 
[13]. 

3. SIMULATION RESULT 

A linear 10 element array is simulated in MATLAB to 
compare results of LMS, Optimized-LMS and RLS 
algorithm. Spacing between adjacent elements of 
array equals one half of the wavelength. 



2.5 
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Figure3. Adaptive beam forming using NLMS 
Algorithm 
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Figure4. Adaptive beam forming using RLS 
Algorithm 
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Figure5. Comparison of SNR vs BER for Adaptive 
beam forming using LMS, NLMS and RLS 
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4. CONCLUSION 

By analyzing the graphs, we observed that RLS 
algorithm provides fastest convergence, Optimized 
LMS algorithm also shows fast convergence but LMS 
algorithm lacks the convergence speed. In beam 
forming results, RLS showed the best beam forming 
capability placing deeper nulls in case of all the three 
interference positions i-e 40±, 60± and 90±. The 
significant difference between the results of LMS and 
Optimized-LMS in case of beam forming was the 
presence of many minor lobes in Optimized-LMS. 
The dependency on SNR and SIR showed that in 
better conditions i.e., high SNR and SIR Optimized- 
LMS showed the best results. But in poor conditions 
i-e low SNR and SIR its performance deteriorates. 
LMS and RLS almost showed equal dependency on 
SNR and SIR. As the recent developments in digital 
signal processor (DSP) kits and field-programmable 
gate arrays (FPGA) have made it possible to 
implement RLS algorithms in real time systems, and 
complexity to an extent is not a problem anymore. So 
RLS algorithm is proposed as it provides deeper 
nulls in the direction of interferences and faster 
convergence. 
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